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Abstract 

Cosmic rays with energies above 10^^ eV are currently of considerable interest 
in astrophysics and are to be further studied in a number of projects which are 
either currently under construction or the subject of well-developed proposals. 
This paper aims to discuss some of the physics of such particles in terms of 
current knowledge and information from particle astrophysics at other energies. 

Keywords: cosmic rays, acceleration of particles, magnetic fields, dark 
matter, radiation mechanisms: nonthermal, diffuse radiation 



1 Introduction 

Cosmic rays (CR) are the non-photon particles with which we learn about astrophysics. 
Their composition ranges over the known nuclei (and antiprotons) to electrons (and 
positrons) and neutrinos of all flavours, and (perhaps) to exotic particles as yet un- 
observed in accelerator physics. Cosmic ray studies are complementary to photon 
astrophysics since many astrophysical photons are produced in processes, such as syn- 
chrotron emission, which involve charged cosmic ray particles. In this review we con- 
centrate on the highest energy cosmic rays, which are unaffected by heliospheric modu- 
lation, but are strongly affected by propagation effects through our galaxy. In the case 
of an extragalacic origin for the highest energy cosmic rays, their propagation through 
intergalactic photon and magnetic fields will have a profound influence on what we 
observe. 

Radio astronomical studies probe cosmic rays at their sources and we are led to 
associate cosmic ray acceleration with energetic radio objects, pulsars, supernovae, 
AGNs etc. The next generation of large radio projects (LOFAR, SKA etc.) will 
certainly add considerably to our understanding of cosmic ray sources, and may even 
offer opportunities to develop new detection techniques through the radio-frequency 
fields of the air showers of cosmic rays arriving at Earth (e.g. Falcke and Gorham 
2003). 



The radio emission is dominated by cosmic ray electrons which rapidly lose energy in 
the emission processes. The more massive nuclei do not have a clear radio signature but 
propagate with little energy loss and are readily observed at Earth. They are believed 
to fill (at least) the volume of our galaxy. However, being charged particles, their paths 
are dictated by astrophysical magnetic fields and only at the highest energies do we 
expect magnetic deflections to be sufliciently small for directional astrophysics to be 
possible. 

At modest energies (up to about lO^^eV), cosmic rays are sufficiently plentiful 
that balloon experiments can have enough collecting area to directly study properties 
of the beam. Indeed, even satellite experiments extend up to 2 TeV/amu (for the 
CRN experiment on Spacelab 2, Miiller et al 1991) which for iron is about 10^'^eV per 
nucleus. At higher energies, we depend on techniques which investigate the cascades of 
particles (extensive air showers) resulting from the impact of cosmic ray particles on our 
atmospheric gas. Since a single cosmic ray particle can produce millions (or billions) 
of secondary particles in this way, and those particles scatter in the atmosphere to 
hundreds of metres from the original cosmic ray trajectory, such studies offer uniquely 
efficient opportunities for studying fluxes down to levels of primary CR particles per 
square kilometre per century. 

The cosmic ray cascades may be directly detected with groups of spaced large-area 
radiation detectors on the ground (ground arrays) or indirectly through collecting (with 
large telescopes) the photons which they produce, fluorescence light from atmospheric 
nitrogen or Cerenkov emission from the bulk atmospheric gas. In the future, such light 
collectors may well be satelhte-based. 

The Ultra High Energy cosmic rays (UHE CR) which arc the focus of this review 
are currently of particular interest. Spectral data at energies above 1 EeV (lO^^eV) 
and directional results, notably from the AC AS A project, are very suggestive of fasci- 
nating, unexpected physics (Hayashida et al. 1999). Furthermore, this fleld of research 
is experimentally challenging and there is controversy in discrepancies between exper- 
imental data from experiments currently operational or recently discontinued. A new 
era in the field is beginning with the commissioning of the Pierre Auger Observatory 
(Auger Collaboration 2001) which comprises a pair of 3000 km^ arrays (one under 
construction in Mendoza Province, Argentina, and one planned for Utah, U.S. A) em- 
ploying both particle and optical detectors, followed by large-scale Japanese projects, 
and possibly the space-based Extreme Universe Space Observatory project (EUSO is 
a Europe/ Japan/US collaboration currently under Phase A study as an ESA mission 
with the goal of a three-year mission on the ISS starting in 2009, see http://www.euso- 
mission.org/). 

At the present time, our ideas concerning cosmic rays at these highest energies are 
predominantly based on the idea that such particles are likely to be of extragalactic 
origin. Their sources are presumed to be found in some extreme astrophysical envi- 
ronment such as the most energetic radio or gamma-ray sources. We have been forced 
to these ideas by our failure to identify any models for galactic objects capable of 
accelerating particles to within one thousandth of the required energies, and also by 
the failure of our galactic source models to reproduce the directional isotropy of the 
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observed beam. Nonetheless, extragalactic scenarios have their own problems. The 
intergalactic magnetic field has largely unknown properties. It could be strong enough, 
with a structure which makes it impenetrable to particles from nearby clusters of galax- 
ies in realistic periods of time. Also, intracluster magnetic fields may limit the ability 
of particles even to leave galactic clusters. 

This paper begins by presenting a brief overview of the observed properties (energy 
spectrum, arrival directions, and composition) of the cosmic ray beam. There follows 
a physical discussion of the cosmic ray acceleration process and resulting limits to 
the maximum achievable energies. Energy losses in the following propagation through 
astrophysical fields arc fundamental in determining the observed beam properties and 
these are then discussed. It is possible that the origin of the highest energy cosmic rays 
hes in exotic particles and that field of research is introduced. Whatever the sources 
of these particles may be, their observed directional properties will be closely linked 
with the properties of the intervening magnetic fields. A discussion of this key factor 
completes our review. For an earlier review see Nagano and Watson (2000). 

1.1 The Detection of Ultra High Energy Cosmic Rays 

Cosmic ray observatories record the arrival directions and distribution of energies of 
incident particles. They also attempt to measure a remaining parameter, which is 
the mass composition of the primary cosmic ray particle. It is in the latter respect 
that cosmic ray studies differ from others in astrophysics. The mass composition is 
a key parameter in our astrophysical understanding since the charge on the particle 
(closely related to mass, for a nucleus), with the momentum (or energy) determines 
the propagation path of the particle. 

At the highest energies, the cosmic ray flux is extremely low and the particle record- 
ing must be through processes which enable detection to be achieved at large distances 
from the path of the primary cosmic ray. In practice, this involves making use of the 
particle cascade which results from the interaction of the primary cosmic ray particle 
with our atmosphere. 

That cascade is initiated by a primary cosmic ray particle when its first atmo- 
spheric interaction occurs. A cascade of secondary particles is then fed by degrading 
energy from the primary particle as it repeatedly interacts in its atmospheric passage. 
Those secondary particles cause energy to be deposited into the atmospheric gas, with 
some remaining energy reaching the groimd. The cosmic ray detection process then 
consists in either a measurement of the passage of the energy as it is carried by parti- 
cles through the atmosphere (through any emitted Cerenkov light or by any induced 
nitrogen fiuorescence light) or in direct detections in radiation detectors of remaining 
particles which reach the ground. 

The arrival direction of the primary particle is deduced from the direction of the 
path of the atmospheric cascade and characteristically has a resolution of a fraction 
of a degree. The primary particle energy is deduced by attempting to make sufficient 
measurements on the cascade such that a calorimetric accounting (to a few tens of 
percent) can be made of the various energy sinks. This is most directly done through 
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nitrogen fluorescence studies, but the signal is weak and the technique works best at the 
highest energies. Mass composition measurements are the most difficult, and cannot 
be made on an individual cascade basis. They depend on statistical studies of the 
cascade developments as a function of energy. Crudely speaking, massive nuclei have a 
large cross section and interact early. Their cascades also degrade in energy relatively 
rapidly. As a result, early developing cascades are signatures of 'heavy' primaries and 
late developing cascades are indicative of 'light' (probably proton) primaries. Massive 
amounts of cascade modelling puts flesh on these arguments. 

Nitrogen fluorescence studies were pioneered by the Fly's Eye experiment (Baltru- 
saitis ct al. 1985) and its successor, the High Resolution Fly's Eye (HiRcs). Many 
ground detector arrays have been used. The largest one in fully developed use is the 
AGASA (Hayashida et al. 1999) experiment in Japan. Both of the HiRes and AGASA 
experiments are sufficiently large to probe energies in excess of lO^^eV. The next of 
these huge experiments will use both techniques. The Pierre Auger Observatory (Auger 
Collaboration 2001) is designed to operate above lO^^cV. It will have a 3000 square kilo- 
metre collecting area instrumented with 1600 large-area particle detectors and will also 
have twenty-four 4 m diameter Schmidt optical telescopes for fluorescence detection. 

1.2 The Cosmic Ray Energy Spectrum 

The cosmic ray energy spectrum is shown in Fig. 1. It is remarkable both in its 
range of energies and in its range of fluxes. It covers over ten decades of energy and 
thirty decades of flux in a form close to a power law with an index of about -2.7. 
Deviations from that power law are relatively small but are generally regarded as 
physically significant. There is a steepening at about lO^^eV, known as the knee, and a 
flattening at about lO^^eV known as the ankle (e.g. Abu-Zayyad et al. 2001). The knee 
is often argued to be associated with an energy limit of acceleration from supernova 
remnant sources, although it may well be related to a loss of ability for our galactic 
magnetic fields to retain (and build up internally) the cosmic ray flux (e.g. Clay 2002). 
The ankle is usually associated with the onset of a dominant, flatter, extra- galactic 
cosmic ray spectrum. It is important to note that, in this model, our galaxy produces 
particles with energies up to those of the ankle of the spectrum. That is already above 
1 EeV=10^^eV, and is well above energies which are easily accessible for any present 
galactic acceleration models. 

A key region of the spectrum is its very highest energies. The flux here is so 
low that the low event statistics in our observations to date leave us uncertain of the 
spectral structure above the key energy of 6 x 10^^ eV. Here there is a predicted spectral 
downturn, the Greisen-Zatsepin-Kuzmin (GZK) cut-off (Greisen 1966, Zatsepin and 
Kuzmin 1966) for particles which have travelled more than a few tens of Mpc, due to 
interactions with the 2.7K cosmic microwave background radiation (CMBR). However, 
several experiments have reported CR events with energies above 10^*^ eV (Takeda et 
al. 2003) with the highest energy event having 300 EeV (Bird et al. 1995). Very recent 
data from the two largest aperture high energy cosmic ray detectors are contradictory: 
AGASA (Takeda et al. 2003) observes no GZK cut-off while HiRes (Abu-Zayyad et 
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Figure 1: The cosmic ray energy spectrum (dotted) as measured from the Earth (after 
Bhattacharjee and Sigl 2000). 
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al. 2002) observes a cut-off consistent with the expected GZK cut-off. A systematic 
over-estimation of energy of about 25% by AGASA or under-estimation of energy of 
about 25% by HiRes could account the discrepancy (Abu-Zayyad et aL 2002), but the 
continuation of the UHECR spectrum to energies well above 10^'' cV is now far from 
certain. Future measurements with Auger (Auger Collaboration 2001) should resolve 
this question. If the spectrum does extend well beyond 10^° eV, determining the origin 
of these particles could have important implications for astrophysics, cosmology and 
particle physics. 

1.3 Arrival Directions 

A key observation in cosmic ray astrophysics is the directional distribution of the 
particles. That distribution will depend on any galactic magnetic fields and hence will 
be energy (rigidity) dependent. However, with very limited exceptions, which arc not 
individually statistically significant, there is no observed deviation from isotropy above 
the knee of the energy spectrum, and any anisotropics at lower energies are themselves 
very small (Smith and Clay 1997, Clay, McDonough and Smith 1997). Recently, the 
AGASA experiment (Hayashida et al. 1999) found a non-uniform distribution of arrival 
directions, suggestive of a source direction, in the energy range lO^^'^^eV to lO'^^'^eV. 
That observation is potentially very important, particularly as there is supporting 
evidence in data from the SUGAR array (Bellido et al. 2001). However, neither of 
those observations on their own is clearly statistically significant. Still, those data are 
regarded by many as the possible beginning of a new era in cosmic ray astrophysics 
in which we can begin directional cosmic ray astronomy. The possibility of having a 
source to observe may indeed open up new frontiers. 

1.4 Composition of the UHECR 

The highest energy cosmic rays show no major differences in their air shower char- 
acteristics to cosmic rays at lower energies. One would therefore expect the highest 
energy cosmic rays to be protons particularly if, as is most likely, they are extragalac- 
tic in origin. However, it is still possible that they are not single nucleons. Obvious 
candidates are heavier nuclei (e.g. Fe), 7-rays and neutrinos. Surprisingly, it is even 
more difficult to propagate nuclei than protons, because of the additional photonuclear 
disintegration (Tkaczyk et al. 1975, Puget et al. 1976, Karakula and Tkaczyk 1993, 
Elbert and Sommers 1995, Anchordoqui et al. 1997, Stecker and Salamon 1999). The 
possibihty that the 300 EeV Fly's Eye event is a 7-ray has been discussed (Halzen et 
al. 1995) and, although not completely ruled out, the air shower development profile 
seems inconsistent with a 7-ray primary. Weakly interacting particles such as neutrinos 
will have no difficulty in propagating over extragalactic distances, of course. This pos- 
sibility has been considered, and generally discounted (Halzen et al. 1995, Elbert and 
Sommers 1995), mainly because of the relative unhkelihood of a neutrino interacting in 
the atmosphere, and the resulting great increase in the luminosity required of cosmic 
sources. 
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Figure 2: Proposed sites for cosmic ray acceleration related to their likely dimensions 
and magnetic field strength for scattering centres moving with Vscat ~ c (after Hillas 
1984). The lines represent plausible limits for 10^° eV cosmic ray containment in the 
sources (solid line - protons, dashed line - iron nuclei) for scattering centres moving 
at at speeds ~ c. For scattering centres moving at speeds Vgcat < c the size would need 
to be a factor c/vscat larger. 

1.5 Cosmic Ray Sources 

Figure 2 is a well known diagram first produced by Hillas (1984). It reminds us that 
acceleration, associated with magnetic structures, requires the field and its dimensions 
to be sufficient to contain the accelerating particle through the acceleration process. 
The lines simply give the magnetic field vs. gyroradii Vg for protons (solid) and iron 
nuclei (dashed), and this gives the minimum size for scattering centres movinff at at 
speeds ~ c. More realistically, for scattering centres moving at speeds f scat < c the size 
would need to be a factor c/fgcat larger. This puts a limit on the product of the source 
field and its physical dimensions. Strong fields with large-scale structure are attractive 
for acceleration to the highest energies. The acceleration is thought to be likely to be 
associated with astrophysical shocks. 

One of the very few plausible acceleration sites may be associated with the radio 
lobes of powerful radio galaxies, either in the hot spots (Rachen and Biermann 1993) or 
possibly the cocoon or jet (Norman et al. 1995). One-shot processes such as magnetic 
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reconnection (e.g. in jets or accretion disks) comprise another possible class of sources 
(Haswell at al. 1992, Sorrell 1987). 

Acceleration at the termination shock of the galactic wind from our Galaxy has also 
been suggested by Jokipii and Morfill (1985), but due to the lack of any statistically 
significant anisotropy associated with the Galaxy it is unlikely to be the explanation. 
However, a re-evaluation of the world data set of cosmic rays has shown that there is a 
correlation of the arrival directions of cosmic rays above 40 EeV with the supergalac- 
tic plane (Stanev et al. 1995), lending support to an extragalactic origin above this 
energy, and in particular to models where "local" sources (< 100 Mpc) would appear 
to cluster near the supergalactic plane. Such a correlation would also be consistent 
with a Gamma Ray Burst (GRB) origin as two type Ic hypernovae (supernovae with 
broad absorption features indicating high velocity ejected material and a rather large 
explosion energy) have now been identified with GRB (SN 2003dh/GRB 030329 Kawa- 
bata et al. 2003; SN 1998bw/GRB980425 Galama et al. 1998). The expanding fireball 
may have ultrarelativistic components (e.g. T ~ 300) and this may lead to production 
of UHECR through relativistic shock acceleration (Vietri 1995) or some other process 
(see e.g. Dermer 2002 for a discussion and references to earher work). 

Because of the resulting fiat spectrum of particles (including gamma-rays and pro- 
tons) extending up to GUT (grand unified theory) scale energies, topological defect 
models have been invoked to try to explain the UHE CR. Propagation of the spectra 
of all particle species over cosmological distances is necessary to predict the cosmic ray 
and gamma-ray spectra expected at Earth. Propagation over cosmological distances 
to Earth (as would be the case in some topological defect origin models) results in 
potentially observable gamma-ray fiuxes at GeV energies in addition to cosmic rays. 
Massive relic particles on the other hand, would cluster in galaxy halos, including that 
of our Galaxy, and may give rise to anisotropic cosmic ray signals at ultra high energies. 

The suggestion that neutron stars might accelerate cosmic rays followed soon af- 
ter the discovery of pulsars (see, e.g. Gold 1975 and references therein). Voltages up 
to ~ 10^^-lO^^V (depending on pulsar period and magnetic field) are expected in a 
pulsar's magnetosphere. These could accelerate nuclei with a resulting fiat spectrum 
extending up to ~10^^ cV, and could possibly explain the knee in the cosmic ray spec- 
trum (e.g. Bednarek & Protheroe 2002). The pulsar wind shock has been proposed as 
an acceleration site (e.g., Berezhko 1994, Bell and Lucek 1996) and might in princi- 
ple accelerate particles to lO^^-lO^^eV. Blasi et al. (2000) suggest that the UHE GR 
events are due to iron nuclei accelerated from young strongly magnetized neutron stars 
through relativistic MHD winds of neutron stars whose initial spin periods arc shorter 
than ~ 10 ms. More recently, the acceleration to ultra-high energies has been discussed 
in the context of fast re-connection in millisecond pulsars, but it appears that to reach 
~ 10^° eV magnetic fields ~ 10^^-10^^ G, and special geometries are required (de Gou- 
veia Dal Pino and Lazarian, 2001). In our opinion, while neutron stars may contribute 
to the observed cosmic ray spectrum it seems unlikely that they arc responsible for the 
UHE CR, with the possible exception of transient shocks in pulsar winds of neutron 
stars formed during a hypernova explosion (i.e. a GRB as already discussed). 

We shall delay a more detailed discussion of radio galaxies, active galactic nuclei 
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and topological defects as sources of the UHE CR, and discuss next diffusive shock 
acceleration and cosmic ray propagation. 



2 Diffusive Shock Acceleration 

For stochastic particle acceleration by electric fields induced by the motion of magnetic 
fields B, the rate of energy gain by relativistic particles of charge Ze can be written 
(in SI units) 



dE 
It 



i{E)Zec'B (1) 



where ^{E) < 1 and depends on the acceleration mechanism. Below is a simple heuristic 
treatment of Fermi acceleration based on those given by Gaisser (1990) and Protheroe 
(2000). We shall start with 2nd order Fermi acceleration (Fermi's original theory) 
and describe how this can be modified in the context of astrophysical shocks into the 
more efficient 1st order Fermi mechanism known as diffusive shock acceleration. More 
detailed and rigorous treatments are given in several review articles (Drury 1983a, 
Blandford and Eichler 1987, Berezhko and Krymsky 1988). See the review by Jones 
and Ellison (1991) on the plasma physics of shock acceleration which also includes a 
brief historical review and refers to early work. 



2.1 Fermi's Original Theory 

Gas clouds in the interstellar medium have random velocities of ~ 15 km/s superim- 
posed on their regular motion around the galaxy. Cosmic rays gain energy on average 
when scattering off these magnetized clouds. A cosmic ray enters a cloud and scatters 
off irregularities in the magnetic field which is tied to the partly ionized cloud. 
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Figure 3: Interaction of cosmic ray of energy Ei with "cloud" moving with speed V 

In the frame of the cloud: (a) there is no change in energy because the scattering is 
coUisionless, and so there is elastic scattering between the cosmic ray and the cloud as 
a whole which is much more massive than the cosmic ray; (b) the cosmic ray's direction 
is randomized by the scattering and it emerges from the cloud in a random direction. 
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Consider an ultra-relativistic cosmic ray entering a cloud with energy Ei and mo- 
mentum pi ^ Ei/c travelling in a direction making angle 9i with the cloud's direction. 
After scattering inside the cloud, it emerges with energy E2 and momentum p2 ~ E2/C 
at angle 62 to the cloud's direction (Fig. 3). The energy change is obtained by applying 
the Lorentz transformations between the laboratory frame (unprimed) and the cloud 
frame (primed). Transforming to the cloud frame: 

E[^-fEi{l- pcosOi) (2) 



where (3 — V/c and 7=1/ Vl — P"^- 
Transforming to the laboratory frame: 

E2^-iE'^{l + (3cose'2). (3) 

The scattering is collisionless, being with the magnetic field. Since the magnetic 
field is tied to the cloud, and the cloud is very massive, in the cloud's rest frame there is 
no change in energy, E2 — E[, and hence we obtain the fractional change in LAB-frame 
energy {E2 -Ei)/Ei, 

AE 1 - COS Oi + 3 cos 0' - /J2 cos 0^ cos 6' 
-^ = ^ ^ 

We need to obtain average values of cos^i and cos ^2- Inside the cloud, the cosmic 
ray scatters off magnetic irregularities many times so that its direction is randomized, 

(cos^^) = 0. (5) 

The average value of cos6'i depends on the rate at which cosmic rays collide with 
clouds at different angles. The rate of collision is proportional to the relative velocity 
between the cloud and the particle so that the probability per unit solid angle of having 
a collision at angle 9i is proportional to {v — l^cos^i). Hence, for ultrarelativistic 
particles {v — c) 

dP 

oc (1 -/3cos^i), (6) 



and we obtain 



givmg 



{AE) 1 + /?V3 ^^4^2 



-Ic^-P' (8) 



E 1-/32 3' 
if /3 < 1. 

We see that {AE) / E oc 13"^ is positive (energy gain), but is 2nd order in (3 and if 
j3 <^ 1 the average energy gain per cloud collision is very small. This is because there 
are almost as many overtaking collisions (energy loss) as there are head-on collisions 
(energy gain). The acceleration rate is 



I dE 



{AE) 4 (V\^ 

-TcoW - o — ^coU (9) 



E ""^^ 3 V c y 
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where Vcow is the rate of colhsion of the cosmic ray with the cloud. 2nd order Fermi 
acceleration is one example of stochastic acceleration. Another example involves the 
average energy gain that arises in resonant interactions of cosmic rays with Alfven 
waves. In this case, the acceleration rate would be 

racc(-E) OC (^^^ TcoU (10) 

where va is the Alfven velocity and Tcou is now the rate of collision of the cosmic ray 
(speed v) with the Alfven waves (see e.g. Jones 1994). 

2.2 1st Order Fermi Acceleration at Astrophysical Shocks 

Fermi's original theory was modified in the 1970's (Axford, Lear and Skadron 1977, 
Krymsky 1977, Bell 1978, Blandford and Ostriker 1978) to describe more efficient 

acceleration (1st order in f3) taking place at supernova shocks but is generally applicable 
to strong shocks in other astrophysical contexts. Our discussion of shock acceleration 
will be of necessity brief, and omit a number of subtleties. 

Here, for simplicity, we adopt the test particle approach (neglecting effects of cosmic 
ray pressure on the shock profile), adopt a plane geometry and consider only non- 
relativistic shocks. Nevertheless, the basic concepts will be described in sufficient detail 
that we can consider acceleration and interactions of the highest energy cosmic rays, 
and to what energies they can be accelerated. We consider the classic example of a SN 
shock. During a supernova explosion several solar masses of material are ejected at a 
speed of ~ 10^ km/s which is much faster than the speed of sound in the interstellar 
medium (ISM) which is ~ 10 km/s. A strong shock wave propagates radially out (speed 
Vs) as the ISM and its associated magnetic field piles up in front of the supernova ejecta 
which moves at speed Vp (see Fig. 4a). As seen in the frame of the shock (see Fig. 4b) 
gas flows from upstream into the shock with speed Ui = Vs and density pi, and flows 
out of the shock downstream with speed U2 — {Vs — Vp) and density p2- The velocity 
of the shock, Vs, depends on the velocity of the piston (ejecta), Vp, and on the ratio of 
specific heats, 7. The compression ratio for a strong non-relativistic shock is given by 

R=Pl = !^ = l±l (11) 

Pl U2 7-1 

from which V5 = Ru2 and Vp = (mi — M2) = {R— 1)^2 giving 

Vs/Vpc^R/{R-l). (12) 

For SN shocks the SN will have ionized the surrounding gas which will therefore be 
monatomic (7 = 5/3), and so a strong shock will have R — A. 

In order to work out the energy gain per shock crossing, we can visualize mag- 
netic irregularities tied to the plasma on either side of the shock as clouds of magne- 
tized plasma of Fermi's original theory (Fig. 5). Here, we assume that the shock is 
non-relativistic, such that we can make the approximation that the ultra-relativistic 
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(a) Upstream frame 



(b) Shock frame 
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Figure 4: A shock moving with speed Vs ahead of shocked gas moving at the pistion 
speed Vp as seen (a) in the upstream frame, (b) in the shock frame. 



accelerated particles are isotropic in both upstream and downstream frames. By con- 
sidering the rate at which cosmic rays cross the shock from downstream to upstream, 
and upstream to downstream, one finds (cos6'i) ~ —2/3 and (cos6'2) ~ 2/3, and Eq. 4 
then gives 

{AE) 4^ AVp A{R-l)Vs 



E ^ 3^ 3 c ~ 3 R c' 
Note this is 1st order in (3 = Vp/c, and so the fractional energy change per collision 
can be much higher than in Fermi's original theory. This is because of the converging 
flow - whichever side of the shock you are on, if you are moving with the plasma, the 
plasma on the other side of the shock is approaching you at speed Vp. 



(13) 




Figure 5: Interaction of cosmic ray of energy Ei with a shock moving with speed Vg. 

To obtain the energy spectrum we need to find the probability of a cosmic ray 
encountering the shock once, twice, three times, etc. If we look at the diffusion of a 
cosmic ray as seen in the rest frame of the shock (Fig. 6), there is clearly a net flow of 
the energetic particle population in the downstream direction. 

The flux of cosmic rays lost downstream is 

/loss = ricnVs/R (14) 
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u^Vs/R 
downstream upstream 

Figure 6: Diffusion of cosmic rays from upstream to downstream seen in tlie sliock 
frame. 

since cosmic rays with number density ncR at the shock are advected downstream 
with speed U2 = Vs/R (from right to left in Fig. 6) and we have neglected relativistic 
transformations of the rates because Vg ^ c. 

Upstream of the shock, cosmic rays travelling at speed v at angle to the shock 
normal (as seen in the laboratory frame) approach the shock with speed (V5 + v cos 9) 
as seen in the shock frame. Clearly, to cross the shock, cos^ > —Vs/v. Then, assuming 
cosmic rays upstream are isotropic, the flux of cosmic rays crossing from upstream to 
downstream is 



1 /■! 

/cross = ^CRi— / (Ks' + ^^cos6l)27rd(cos6') fti ncRv/4:. (15) 
47r J -Vs/v 



The probability of crossing the shock once and then escaping from the shock (being 
lost downstream) is the ratio of these two fluxes: 

Prob. (escape) = /ioss//cross ~ Ws/Rv. (16) 

The probability of returning to the shock after crossing from upstream to downstream 
is 

Prob. (return) — 1 — Prob. (escape), (17) 

and so the probability of returning to the shock m times and also of crossing the shock 
at least m times is 

Prob. (cross > m) = [1 — Prob. (escape)]*". (18) 
The energy after m shock crossings is 

E^Eoi^l + j (19) 

where Eq is the initial energy. 

To derive the spectrum, we note that the integral energy spectrum (number of 
particles with energy greater than E) on acceleration must be 

N{> E)<x[l- Prob.(escape)]'" (20) 
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where 

m - (21) 
ln{l + AE/Ey 



Hence, 



ln7V(> E)^A+ ^^^ff^^j^^ ln[l - Prob. (escape)], (22) 



where ^4 is a constant, and so 

\nN{> E) ^ B - (T -l)\nE (23) 

where i? is a constant and 

/ _ ln[l-Prob.(escape)] \ _ R + 2 

' ln(l + AE/E) J ^ R-1 ^ ^ 



where we have used ln(l + x) ~ x for a; <^ 1. 

Hence we arrive at the spectrum of cosmic rays on acceleration 

A^(> E) oc (integral form) (25) 

^ oc E-^ (differential form). (26) 
dE 

For compression ratio R — A (strong shock) we have the well-known E"'^ differential 
spectrum. The observed cosmic ray spectrum is steepened by energy-dependent escape 
of cosmic rays from the Galaxy. 

2.3 Acceleration Rate 

Here we again neglect effects of cosmic ray pressure and consider only a non-relativistic 
shock. The acceleration rate is defined by 



.(E) . i 
^ E dt 



{{^E)/E) ^ 4 (it:-l) Vs ^27) 

^cycle(-E') 3 R Ctcycle{E) 



where tcycie is the time for one complete cycle, i.e. from crossing the shock from up- 
stream to downstream, diffusing back toward the shock and crossing from downstream 
to upstream, and finally returning to the shock. 

The rate of loss of accelerated particles downstream is the probability of escape per 
shock crossing divided by the cycle time 

, Prob. (escape) A Vs 3 , , 

re.o{E) = ^ ^ ^ - ' = -—r^.{E) (28) 

''cycle -ft C t-cycle ^.-'-'^ 

We see immediately that the ratio of the escape rate to the acceleration rate depends 
on the compression ratio 
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and for a strong shock (R = 4) the two rates are equal, giving the well-known E~'^ 
power-law. 

We shall discuss these processes in the shock frame (see Fig. 6) and consider first 
particles crossing the shock from upstream to downstream and diffusing back to the 
shock, i.e. we shall work out the average time spent downstream. Since we are con- 
sidering non-relativistic shocks, the time scales are approximately the same whether 
measured in the upstream or downstream plasma frame, and so in this section we drop 
the use of subscripts indicating the frame of reference. 

Diffusion takes place in the presence of advection at speed U2 in the downstream 
direction. The diffusion coefficient will be a function of magnetic rigidity which, for 
ultra-relativistic particles considered in this paper, is approximately equal to E/Ze 
where Ze is the charge. However, here we are mainly concerned with singly charged 
particles and shall work in terms of E rather than rigidity. The diffusion coefficient 
depends on the turbulence in the magnetic field. Often a power-law dependence k{E) oc 
E^ is assumed, where S depends on the spectrum of turbulence. For a Kolmogorov 
spectrum 6 = 1/3, and for a completely tangled magnetic field S = 1 at least over some 
range of energies. 

The typical distance a particle diffuses in time t is -\/^ where k2 is the diffusion 

coefficient in the downstream region. The distance advectcd in this time is simply U2t. 
If \/k2t ^ U2t the particle has a very high probability of returning to the shock, and 
if \/k2t -C U2t the particle has a very high probability of never returning to the shock 
(i.e. it has effectively escaped downstream). So, we set ^/k2t = U2t to define a distance 
k2/u2 downstream of the shock which is effectively a boundary between the region 
closer to the shock where the particles will usually return to the shock and the region 
farther from the shock in which the particles will usually be advected downstream 
never to return. There are ncRk2/u2 particles per unit area of shock between the shock 
and this boundary. Dividing this by /cross (Eq. 15) we obtain the average time spent 
downstream before returning to the shock 

^ ifeW. (30) 

C U2 

Consider next the other half of the cycle after the particle has crossed the shock 
from downstream to upstream until it returns to the shock. In this case we can define 
a boundary at a distance ki/ui upstream of the shock such that nearly all particles up- 
stream of this boundary have never encountered the shock, and nearly all the particles 
between this boundary and the shock have diffused there from the shock. Then dividing 
the number of particles per unit area of shock between the shock and this boundary, 
ncRki/ui, by /cross we obtain the average time spent upstream before returning to the 
shock 

4fci(E) 



and hence the cycle time 



ti{E) ^ (31) 

C Ui 



1(M^ + M^). (32) 

c V Ml U2 I 
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The acceleration rate is then given by 




-1 



racc(-E') 



(33) 



At this point, a comparison with stochastic acceleration is appropriate. Noting that 
the diffusion coefficient can be written 



where Acou and Tcou are the effective collision mean free path and collision frequency, 
respectively, the acceleration rate can be written 



which has the same functional form as for stochastic acceleration acceleration (Eq. 10) 
as pointed out by Jones (1994). 

For either shock acceleration or stochastic acceleration to be able to accelerate 
cosmic rays to high energies the physical conditions must be suitable: Alfven waves, 
magnetic irregularities or turbulence in the magnetic field must be present on length 
scales of the gyroradii of particles being accelerated and provide a sufficiently high 
scattering rate such that the required maximum energy can be achieved during the life- 
time of the accelerator, and the large-scale magnetic field must be able to confine the 
the highest energy particles within the accelerator. There are only a few places (solar 
fiares) where the turbulence is likely to be strong enough for stochastic acceleration to 
work, and the spectrum will differ for different species. Generally, shock acceleration 
is favoured for the following reasons: (a) coUisionless shocks exist everywhere, provide 
the necessary physical conditions, and are known to accelerate particles efficiently; (b) 
the energy associated with shocks can be large (a significant fraction of the energy 
released in a supernova explosion is carried by the supernova ejecta); (c) in the test 
particle case at least, the spectrum of accelerated particles is a power law which is the 
same for all species and depends only on the compression ratio (see Eq. 24). 

2.4 Maximum Acceleration Rate 

Wc next consider the diffusion for the cases of parallel, oblique, and perpendicular 
shocks, and estimate the maximum acceleration rate for these cases. The diffusion 
coefficients required, ki and k2, are the coefficients for diffusion parallel to the shock 
normal. The diffusion coefficient along the magnetic field direction is some factor 77 
times the minimum diffusion coefficient, known as the Bohm diffusion coefficient. 




(34) 




(35) 



^1 




(36) 



where Tg is the gyroradius, and 77 > 1. 
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Parallel shocks are defined such that the shock normal is parallel to the magnetic 
field (i^lliTi). In this case, making the approximation that ki — k2 = k\\ and Bi = B2 
one obtains 

20 r]E 
3 eBiu{ 

For a shock speed of Ui = 0.1c and 77 = 10 one obtains an acceleration rate (in SI units) 
of 

dE 



II ' /on\ 



dt 



1.5 X IQ-^ec^B. (3^ 



For the oblique case, the angle between the magnetic field direction and the shock 
normal is different in the upstream and downstream regions, and the direction of the 
plasma flow also changes at the shock. The diffusion coefficient in the direction at 
angle 6 to the magnetic field direction is given by 

k = A;|| cos^ e + k± sin^ 6 (39) 

where k± is the diffusion coefficient perpendicular to the magnetic field. Jokipii (1987) 
shows that 

k± (40) 

provided that rj is not too large {rj values up to ~10), and that acceleration at perpen- 
dicular shocks can be much faster than for the parallel case. For a perpendicular shock 
{9 — 90°) , k — k± and B2 ~ 4Si and one obtains 

^^^"^l^i^,- ^^^^ 

For a shock speed of ui — 0.1c and 77 = 10 one obtains an acceleration rate (in SI units) 
of 



dE 
'dt 



OMec^B. (42) 



EUison, Baring and Jones (1995) have examined in detail the acceleration time 
and injection efficiency for oblique shocks. They point out that Eq. 33 is only valid 
when ki/ui > Vg which requires that cosmic ray velocities, v, satisfy r] < v/ui. If this 
condition is not met, as would typically be the case for thermal particles ifrj ^ 1, then 
these particles would have a reduced probability of returning to the shock after crossing. 
As injection into the acceleration process is assumed to be from the thermal plasma, 
having particle speeds v ^ ui, this will cause problems for the injection of cosmic rays 
in the case of highly oblique shocks. Thus, although having 77 3> 1 can significantly 
reduce the acceleration time in oblique shocks, the injection efficiency would also be 
significantly reduced. Ellison, Baring and Jones (1995) find this to be a serious problem 
for 9 > 40°. One possible solution to this problem could be injection of a previously 
accelerated particle population, e.g. acceleration in a pulsar magnetosphere followed 
by injection into the supernova shock. Alternatively, injection from a supra-thermal 
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tail of the thermal distribution could help, with the supra-thermal tail being due to 
heating by hard X-rays or gamma-rays from a nearby source. 

Supernova shocks remain strong enough to continue accelerating cosmic rays for 
about 1000 years. The rate at which cosmic rays are accelerated is inversely propor- 
tional to the diffusion coefficient (faster diffusion means less time near the shock). For 
the maximum feasible acceleration rate, a typical interstellar magnetic field, and 1000 
years for acceleration, energies of 10^^ x Z eV are in principle possible {Z is atomic 
number) at parallel shocks, and 10^^ x Z eV at perpendicular shocks but in this case, 
as noted above, injection may be a problem. 

2.5 Effect of Cosmic Ray Pressure 

Inclusion of the effects of cosmic ray pressure on the shock profile, and consequently 
on the spectrum of accelerated particles, is a very difficult problem and we refer the 
reader to Jones & Ellison (1991) for a detailed discussion. The shock profile, instead of 
being a step-function becomes smoothed, and this affects the acceleration of low and 
high energy particles differently thereby affecting the cosmic ray spectrum. Generally, 
the results are sensitive to Mach number, fraction of energy flux of upstream plasma 
converted to accelerated particles, energy dependence of diffusion coefficient, etc. Cal- 
culated spectra may be E~'^, or flatter, or concave, depending on the input and the 
approximations made. 

The original method of treating this non-linear effect is the two-fluid model (see e.g. 
the review by Drury 1983b), the two fluids being plasma with ratio of speciflc heats 
7=5/3, and cosmic rays with 7=4/3. To solve the steady-state two-fluid equations it is 
necessary to make some approximations, usually that the cosmic rays interact with the 
gas only through their pressure, and that the cosmic ray pressure and energy flux are 
continuous across the shock. Furthermore an effective ratio of speciflc heats 7eff and 
energy-independent diffusion was generally assumed. The solutions were often found to 
be unstable at high Mach numbers unless the spectrum was cut-off artificially, and the 
approximations used meant that there was, in effect, injection without conservation of 
particle number. 

In time-dependent two-fluid model calculations (e.g. Falle & Giddings 1987, Bell 
1987), 7efr can be calculated self-consistently by weighting 7 by the pressures of the 
two fluids taking account of the spectrum of energetic particles, and energy- dependent 
diffusion can be included. Also, finite times for acceleration effectively eliminate the 
problem of injection without conservation of particle number. See Baring (1997) for 
a brief review and additional references, and Blandford and Eichler (1987), Berezhko 
(1999) and Berezhko and Volk (2000) for alternative approaches to this non-linear 
problem. 

2.6 Relativist ic Shocks 

Shocks in jets of active galactic nuclei (AGN) and GRB are hkely to be relativistic, 
i.e. ui > 0.1c, with shock Lorentz factors, F ~ 10 and F ~ 300, respectively. For 
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relativistic plasma motion with bulk velocities comparable to those of the particles 
being accelerated, the particle distribution will not be isotropic, and the approximations 
used earlier are no longer valid. Instead, the typical escape probability and fractional 
energy gain per shock crossing are (very crudely) ~ 0.5 and ~ 1, respectively. Initial 
work on relativistic shock acceleration was done by Peacock (1981); see Kirk & Duffy 
(1999) for a topical review and additional references. 
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Figure 7: Momentum spectra obtained from Monte Carlo simulations with small angle 
scattering upstream for Fi = 10 and k (x. p [5 = 1). Spectra of particles which have 
escaped at times t' /t^^ — 10°, 10^, . . . , 10^^ after injection are shown (sohd histograms) 
together with the spectra of particles remaining in the acceleration zone (dotted his- 
tograms) at these times (adapted from Protheroe 2001). 

The techniques used have been analytic (eigenfunction method) and Monte Carlo 
methods to solve the steady state equations for the particle spectral and angular distri- 
butions. The acceleration time depends strongly on Ui, k±/k\\, and 9 and can be as low 
as ~ (1 — 10)rg/c (e.g. Bednarz & Ostrowski 1996, 1998, and Ostrowski 1998). Detailed 
studies have shown a trend in which shocks with larger ui generally have lower spec- 
tral indices (Kirk and Schneider 1987, Ellison and Jones 1990), and the spectral index 
can be very sensitive to the pitch angle particularly for mildly-relativistic shocks (see 
Baring 1997 for additional references). Nevertheless, for plane ultra-relativistic shocks, 
test-particle Monte Carlo simulations tend to give spectral indices close to F = 2.25 
(e.g. Bednarz & Ostrowski 1998, Gallant, Achterberg & Kirk 1999, Bednarz 2000, Kirk 
et al. 2000, Achterberg ct al. 2001, Protheroe 2001, Meh and Quenby 2003a,2003b). 
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These simulations assume strong turbulence downstream of the shock, but if this not 
present then the spectral index will be steeper (Ostrowski & Bednarz 2002). 

Shock modification by the back-reaction of accelerated particles (EUison & Double 
2002) can cause the compression ratio to increase above the test particle value causing 
the spectrum of accelerated particles to differ from a simple power law with F !=a 2.25 
for shock Lorentz factors less than 10. 

A recent development has been to simulate separately the propagation upstream 
and downstream to work out the probability of returning to the shock at a particular 
angle to the shock normal for a given direction on crossing the shock (Protheroe 2001, 
Lemoine & Pelletier 2003), and use these distributions to simulate very efficiently 
relativistic shock acceleration over a large dynamic range in particle momentum. The 
time evolution of the momentum spectrum for injection at the shock, at time t — 0, of 
highly relativistic mono-energetic particles which are isotropic in the upstream frame 
and have upstream-frame momentum pq as calculated by Protheroe (2001) is shown 
in Fig. 7. Here, as observed in the downstream frame these injected particles have a 
range of initial momenta distributed up to ~ 2pq where Pq = F12P0 (ri2 is the Lorentz 
factor for transforming between upstream and downstream frames) . All distances are 
measured in units of the upstream scattering mean free path and all times are 
measured in units of = X^/c. The figure shows how the spectrum develops with 
time after injection, showing separately at each time indicated the spectrum of particles 
which have escaped, and (dotted curves) the spectrum of particles remaining within 
the acceleration zone (i.e. not having yet escaped downstream). 

3 Interactions of High Energy Cosmic Rays 

Interactions of cosmic rays with radiation are important both during acceleration when 

the resulting energy losses compete with energy gains by, for example, shock acceler- 
ation, and during propagation from the acceleration region to the observer. For UHE 
CR the most important processes are pion photoproduction and Bethe-Heitler pair 
production both on the CMBR, and synchrotron radiation. In the case of nuclei, 
photodisintegration on the CMBR is also important. 

3.1 Nucleons 

The mean interaction length, Xp^, of a proton of energy E is given by, 

1 1 r°° n(e) /■«max(£,-B) „ ^ 

T1fS = ?SF^ r..^/ 'y{s){s-my)dsde^ (43) 

where n{e) is the differential photon number density of photons of energy e, and <j{s) is 
the appropriate total cross section for the process in question for a centre of momentum 
(CM) frame energy squared, s, given by 

s^my + 2sE(l-(3cose) (44) 
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where is the angle between the directions of the proton and photon, and (5c is the 
proton's velocity. For pion photoproduction Smin ~ 1-16 GeV^, and for Bethe-Heitler 
pair production the threshold is somewhat lower, Smin ~ 0.882 GeV^. For both pro- 
cesses, Smax ~ ("^pC^ + AeE) which corresponds to a head-on coUision of a proton of 
energy E and a photon of energy e. 

The mean interaction lengths for both processes, Xjrf{E), are obtained from Equa- 
tion 43 for interactions in the CMBR and are plotted as dashed hues in Fig. 8(a). 
Dividing by the mean inelasticity of the coUision, k,{E), one obtains the energy- loss 
distances for the two processes (sohd curves) , 



E _ Xp^{E) 



dE/dx n{E) 



(45) 
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Figure 8: (a) Mean interaction length (dashed lines) and energy-loss distance (solid 
lines), E/ (dE/dx), for Bethe-Heitler pair production and pion photoproduction in the 
CMBR (lower and higher energy curves respectively). (From Protheroe and Johnson 
1995). (b) Energy-loss distance of Fe- nuclei in the CMBR for pair- production (leftmost 
dashed line) and pion photoproduction (rightmost dashed line). Photodisintegration 
distances are given for loss of one nucleon (lower dotted curve), two nuclcons (upper 
dotted line) as well as the total loss distance (thick curve) estimated by Stcckcr and 
Salamon (1999). The thin full curve shows an estimate over a larger range of energy 
(Protheroe, unpublished) of the total loss distance based on photodisintegration cross 
sections of Karakula and Tkaczyk (1993). 
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3.2 Nuclei 



In the case of nuclei the situation is a httle more comphcated. The threshold condition 
for Bethe-Heitler pair production can be expressed as 




(46) 

and the threshold condition for pion photoproduction can be expressed as 

Since 7 = E/AmpC^, where A is the mass number, we will need to shift both energy-loss 
distance curves in Fig. 8(a) to higher energies by a factor of A. We shall also need to 
shift the curves up or down as discussed below. 

For Bethe-Heitler pair production the energy lost by a nucleus in each coUision near 
threshold is approximately AE f» ^2me(?. Hence the inelasticity is 

2me , 
E Arup ^ ' 

and is a factor of A lower than for protons. On the other hand, the cross section 
goes hke Z^, so the overall shift is down (to lower energy-loss distance) by Z"^ /A. For 
example, for iron nuclei the energy loss distance for pair production is reduced by a 
factor 26V56 12.1. 

For pion production the energy lost by a nucleus in each collision near threshold is 
approximately /S.E a; ■yrriT^c^, and so, as for pair production, the inelasticity is factor 
A lower than for protons. The cross section increases approximately as A^-^ giving an 
overall increase in the energy loss distance for pion production of a factor A^-^ 1.5 
for iron nuclei. The energy loss distances for pair production and pion photoproduction 
are shown for iron nuclei in Fig. 8(b) 

Photodisintegration can be very important both during acceleration and propaga- 
tion and has been considered in detail by Tkaczyk et al. (1975), Puget et al. (1976), 
Karakula and Tkaczyk (1993), Epele and Roulet (1998) and Stecker and Salamon 
(1999). The photodisintegration distance defined by A/(dA/dx) taken from Stecker 
and Salamon (1999) is shown in Fig. 8(b) together with an estimate made over a 
larger range of energy by Protheroe (unpublished) of the total loss distance based on 
photodisintegration cross sections of Karakula and Tkaczyk (1993). Clearly, photodis- 
integration is the dominant loss process for iron nuclei. 



4 Maximum Energies 

Because of their much lower energy losses at a given energy, protons and nuclei can be 
accelerated to much higher energies than electrons for a given magnetic environment. 
For stochastic particle acceleration by electric fields induced by motion of magnetic 



22 



fields B, the rate of energy gain by relativistic particles of charge Ze can be written (in 
SI units) dE/dt = ^Zec^B as in Eq. 1, where ^ < 1 and depends on the acceleration 
mechanism; a value of ^ = 0.04 might be achieved by first order Fermi acceleration at 
a perpendicular shock with shock speed ~ 0.1c. 

The rate of energy loss by synchrotron radiation of a particle of mass Anip, charge 
Ze, and energy ■ymc'^ is 

dE 4 f Z^nieV B^ 2 

7- = -(^T — - 7 c. (49) 

dt 3 ' V y 2//0 ' ^ ' 

Equating the rate of energy gain with the rate of energy loss by synchrotron radiation 
places one limit on the maximum energy achievable by electrons, protons and nuclei: 

GeV, (50) 
GeV, (51) 

1/2 

GeV. (52) 

The cut-off energies of protons and iron nuclei allowed by synchrotron radiation losses 
are shown in Figs. 9(a) and 9(b) respectively, and are plotted against magnetic field 
for three values of ^. 

Equating the total energy loss rate for proton-photon interactions (i.e. the sum of 
pion production and Bcthc-Hcitler pair production) in Fig. 8(a) to the rate of energy 
gain by acceleration gives the maximum proton energy in the absence of other loss 
processes. This is shown in Fig. 9(a) for the three ^ values. As can be seen, for a 
perpendicular shock it is possible in principle to accelerate protons to ~ 10^^ GeV in 
a ~ 10-5 G field. 

The effective loss distance given in Fig. 8(b) is used together with the acceleration 
rate for iron nuclei to obtain the maximum energy as a function of magnetic field. This 
is shown in Fig. 9(b). We see that for a perpendicular shock it is in principle possible 
to accelerate iron nuclei to ~ 2 x 10^^ GeV in a ~ 3 X 10 ^ G field. While this is 
higher than for protons, iron nuclei are likely to get photodisintegrated into nucleons 
of maximum energy ~ 4 x 10^^ GeV, and so there is not much to be gained unless the 
source is nearby. 

Of course, potential acceleration sites need to have the appropriate combination of 
size (much larger than the gyroradius at the maximum energy), magnetic field, shock 
velocity (or other relevant velocity) and the time available for acceleration. These limits 
were obtained and discussed in some detail by Biermann and Strittmattcr (1987). We 
emphasize that the cut-off energies estimated above apply to ideal conditions in which 
strong scattering occurs at all energies between injection and the cut-off energy; in 
reality, the ^-values above are probably optimistic. 



= 6.0x10^/^1^^^ 

. n s -1/2 

= 2.0 X lo'^e^/^ (^) 
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Figure 9: Maximum energy as a function of magnetic field of (a) protons and (b) iron 
nuclei for maximum possible acceleration rate ^ = 1 (dashed), plausible acceleration 
at perpendicular shock ^ = 0.04 (solid), and plausible acceleration at parallel shock 
^ = 1.5 X 10~^ (dotted). Straight hnes on the right give the limit from synchrotron 
loss, curved lines on the left give the limit from Bethe-Heitler pair production and pion 
photoproduction (protons) or photodisintegration (iron nuclei). 



5 Spectral Shape near Maximum Energy 

To determine the spectral shape near maximum energy we use the leaky-box acceler- 
ation model (Szabo and Protheroe 1994, Protheroe and Stanev 1998) which may be 
considered as follows. A particle of energy Eq is injected into the leaky box. While 
inside the box, the particle's energy changes at a rate dE/dt = Era,cc{E) and in any 
short time interval At the particle has a probability of escaping from the box given by 
AtresciE). The energy spectrum of particles escaping from the box then approximates 
the spectrum of shock accelerated particles. 

At time t after injecting Nq particles of momentum Eq into the acceleration zone, 
the number of particles remaining in the acceleration zone, N,.g^{t), is obtained by 
solving 



dt 



which has solution 



Nrem{t) = Nq exp 



- fr^,[E{t)]dt 
Jo 



NqBxp 



The spectrum of accelerated (escaping) particles is then 



dN 
dE 



dN,,^ dt N,,^[t{E)]r^,[E{t)] 



dt dE 



Er,,,[E{t)] 
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rcscjE) 



dE 



(53) 



Let us consider first the case of no energy losses, interactions, or losses due to any 
other process. Assuming that the difi^usion coefficients upstream and downstream have 
the same power-law dependence on energy, and using Eq. 29, 



race = aE-\ resc = (r - l)aE-^ . 



(54) 



Then the differential energy spectrum of particles which have escaped from the accel- 
erator is given by 



dN/dE = No{T-l){Eo)-\E/Eo)-^, 
for {E > Eo) where F — {R + 2)/{R — 1) is the differential spectral index. 



(55) 



5.1 Cut-Off due to Finite Acceleration Volume, etc. 

Even in the absence of energy losses, acceleration usually ceases at some energy due to 
the finite size of the acceleration volume (e.g. when the gyroradius becomes comparable 
to the characteristic size of the shock), or as a result of some other process. We 
approximate the effect of this by introducing a constant term to the expression for the 
escape rate: 

rese = (r - l)aE-' + (r - l)aE-i,. (56) 

where ii^max is defined by the above equation and will be close to the energy at which 
the spectrum steepens due to the constant escape term. We shall refer to E^ax as the 
"maximum energy" even though some particles will be accelerated to energies above 
this. 

Following the same procedure as for the case of a purely power-law dependence of 
the escape rate, one obtains the differential energy spectrum of particles {E > Eq) 
escaping from the accelerator, 

^ = No{r-l){Eo)-\E/Eo)-'[l + {E/E^^)'] 

for 6 > (Protheroe and Stanev 1998). We compare in Fig 10(a) the spectra for F = 2 
and 5 ranging from 1/3 to 1, and note that the energy dependence of the diffusion 
coefficient has a profound influence on the shape of the cut-off. This smooth cut-off 
occurs over up to three decades in energy and its shape depends on the momentum 
dependence of the diffusion coefficients, and the intrinsic spectral index F which de- 
pends on the compression ratio. Such a smooth cut-off has very recently been noted in 
test-particle Monte Carlo simulations of shock acceleration at shocks in a cylindrical jet 
geometry where there is sideways leakage out of the jet (Casse and Markowith, 2003). 
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Figure 10: (a) Differential energy spectrum for tlie case of a cut-off due to escape for 
r = 2 and 6 = 1/3 (solid curve), 2/3 (dotted curve) and 1 (dashed curve) - from 
Protheroe & Stanev (1998). (b) Differential energy spectrum for the case of a cut-off 
due to E"^ energy losses for F = 2 and 5 — 1/3 and -E'cutZ-E'max = 10~^ (leftmost curve), 
10~^, . . . , 10^ (rightmost curve). 



5.2 Cut-Off due to Energy Losses 

When continuous energy losses are included the spectrum is cut-off sharply at an energy 
at which the total rate of energy gain is zero. Depending on the spectral index, and 
momentum dependence of the diffusion coefficient, either a pile- up or a steepening 
in the spectrum occurs before the cut-off. To calculate the energy spectrum a term 
representing the energy-loss rate must be added to the acceleration rate. 



race ^aE ^ + — 



dE 
It 



(58) 



but since the physical size of the "box" increases with energy, synchrotron losses can 
cause a particle in the downstream region to effectively fall out of the box (Drury et al 
1999). This process can be represented by an additional escape term in the escape rate. 
The acceleration zone extends distances Li{E) = ki{E)/ui and L2{E) = k2{E)/u2 
upstream and downstream of the shock. With the additional term in the rate of escape 
of particles due to energy loss. 



(59) 



resc = {T-l)a{E-' + E-i,) + 5-^n,,,{E) (60) 
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where rioss(-E') = —{dE/dt)ioss/E, and we adopt L2/{Li + L2) = R/{1 + R) following 
Drury et al. (1999). 

For the case of synchrotron losses the result depends on the parameters 5, F, Eq, 
Ecut and £^max- As a result of the energy loss by particles near the cut-off energy, a 
pile-up in the spectrum may be produced just below E^ut- The size of the pile-up will 
be determined by the relative importance of race and r^sc at energies just below -Ecut- 
Numerical solution of Eq. 53 gives the results for F = 2, 5 = 1/3 and various -Emax 
which are shown in Fig. 10(b). 

One can use the Monte Carlo method to investigate the shape of the cut-off or pile- 
up which results when the nominal cut-off energy is determined by interactions rather 
than continuous energy losses. This technique was used by Protheroe and Stanev 
(1998) to investigate cut-offs in electron spectra due to inverse Compton scattering 
in the Klein-Nishina regime, and by Szabo and Protheroe(1994) to investigate cut- 
offs in the proton spectrum due to photoproduction in a radiation field. Results for 
E^^t = 2 X 10^2 (j^g photoproduction on the CMBR are shown in Figure 11. 
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Figure 11: The spectrum of particles produced during acceleration (multiplied by E'^) 
per proton injected into the accelerator (T = 2, k oc E): protons (full curve), neutrons 
(dotted curve), charged pions (dashed curve) and neutral pions (chain curve). Results 
arc shown for E^ut = 2 x 10^^ GeV due to photoproduction on the cosmic microwave 
background. (Adapted from figure 7 of Szabo and Protheroe 1994). 



6 Cascading in Cosmic Radiation Fields 

As well as particles being produced during the acceleration process as a result of inter- 
actions, during propagation to Earth cascading occurs and the accompanying fluxes of 
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7-rays and neutrinos must not exceed the observed flux or flux limits. By measuring 
the accompanying fluxes, we may well provide additional clues to the nature and origin 
of the highest energy cosmic rays (Waxman and Bahcall 1999, Mannheim et al. 2001). 
Hence it is important to calculate these fluxes resulting from cascading. 
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Figure 12: (a) The mean interaction length for pair production for 7-rays in: radio 
background calculated by Protheroc and Biermann 1996 (solid curves labelled R), 
radio background of Clark (1970) (dotted line); CMBR (2.7K), infrared and optical 
background (IR) (Malkan and Stecker 1998). The mean interaction length for double 
pair production (4e) in the CMBR (Protheroc and Johnson 1995) are also shown, 
(b) The mean interaction length (dashed line) and energy-loss distance (solid line), 
E/{dE/dx), calculated by Protheroc and Johnson (1995) for electron-photon triplet 
pair production (TPP) and inverse- Compton scattering (IC) in the CMBR. The energy- 
loss distance for synchrotron radiation is also shown (dotted lines) for intergalactic 
magnetic fields of 10"^ (bottom), 10"^°, 10-^\ and 10"^^ gauss (top). 

There are several cascade processes which are important for UHE CR propagating 
over large distances through a radiation field: protons interact with photons resulting 
in pion production and pair production; electrons interact via inverse- Compton scat- 
tering and triplet pair production, and emit synchrotron radiation in the intergalactic 
magnetic field; 7-rays interact by photon-photon pair production. Energy losses due 
to cosmological redshifting of high energy particles and 7-rays can also be important, 
and the cosmological redshifting of the background radiation fields means that energy 
thresholds and interaction lengths for the above processes also change with epoch (see 
e.g. Protheroc et al. 1995). 

The energy density of the extragalactic background radiation is dominated by the 
CMBR. Other components of the extragalactic background radiation are discussed in 
the review of Ressel and Turner (1990). The extragalactic radiation fields which are 
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important for cascades initiated by UHE cosmic rays include the cosmic microwave 
background, the radio background and the infrared-optical background. The radio 
background was measured over thirty years ago (Bridle 1967, Clarke et al. 1970), but 
the fraction of this radio background which is truly extragalactic, and not contamina- 
tion from our own Galaxy, is still debatable. Berezinsky (1969) was first to calculate 
the mean free path on the radio background. More recently Protheroe and Biermann 
(1996) have made a new calculation of the extragalactic radio background radiation 
down to kHz frequencies. The main contribution to the background is from normal 
galaxies and is uncertain due to uncertainties in their evolution. The mean free path 
of photons in this radiation field as well as in the microwave and infrared backgrounds 
is shown in Fig. 12(a). 

Inverse Compton interactions of high energy electrons and triplet pair production 
can be modelled by the Monte Carlo technique (e.g. Protheroe 1986, Protheroe 1990, 
Protheroe ct al. 1992, Mastichiadis et al. 1994), and the mean interaction lengths and 
energy-loss distances for these processes are given in Fig. 12(b). Synchrotron losses 
must also be included in calculations and the energy-loss distance has been added to 
Fig. 12(b) for various magnetic fields. 

Where possible, to take account of the exact energy dependences of cross-sections, 
one can use the Monte Carlo method. However, direct application of Monte Carlo 
techniques to cascades dominated by the physical processes described above over cos- 
mological distances takes excessive computing time. Another approach based on the 
matrix multiplication method has been described by Protheroe (1986) and developed 
in later papers (Protheroe and Stanev 1993, Protheroe and Johnson 1995). A Monte 
Carlo program is used to calculate the yields of secondary particles due to interactions 
with radiation, and spectra of produced pions are decayed to give yields of 7-rays, elec- 
trons and neutrinos. For the pion photoproduction interactions a new program called 
SOPHIA is available (Miicke ct al. 1998). 

7 Radio Galaxies and Active Galactic Nuclei 

Rachen and Biermann (1993) have demonstrated that cosmic ray acceleration hotspots 
of giant radio lobes of Fanaroff- Riley Class II radio galaxies can fit the observed spectral 
shape and the normalization at 10 - 100 EcV to within a factor of less than 10. 
Protheroe and Johnson (1995) repeated Rachen and Biermann's calculation to calculate 
the flux of diffuse neutrinos and 7-rays which would accompany the UHE cosmic rays 
as a result of pion photoproduction on the CMBR, and their calculated flux is shown 
in Fig. 13. The fiux of extremely high energy neutrinos may give important clues to 
the origin of the UHE cosmic rays (for reviews of high energy neutrino astrophysics see 
Protheroe 1998 and Learned and Mannheim 2000). They may even be able to produce 
the observed UHE CR above the GZK threshold through interacting with cosmological 
neutrinos in our galactic halo as discussed in the next section on "Z-bursts" . 

AGN jets may also accelerate protons to ultra high energies and produce neutrino, 
gamma-ray and cosmic ray signals as a result of pion photoproduction interactions in 
the intense AGN radiation fields. There are different versions of these models in which 
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Figure 13: Cosmic ray proton intensity multiplied by E"^ in the model of Rachen and 
Biermann (1993) as calculated by Protheroe and Johnson (1995) for proton injection 
up to 3 X 10^^ GeV (solid line). Also shown are intensities of neutrinos (dotted lines, 
z/^, i/^, z/e, i>e from top to bottom), and photons (long dashed lines). Data are from 
Gaisser and Stanev (1998). 



the target photons are produced inside the blob, e.g. as synchrotron emission by a co- 
accelerated population of electrons (Mannheim 1993, 1995), or are external to the jet, 
e.g. from an accretion disk (Protheroe 1997). In addition, proton synchrotron blazar 
models in which the high energy part of the spectral energy distribution is mainly due 
to synchrotron radiation by protons have been proposed for some blazars (Miicke and 
Protheroe 2000, Miicke et al. 2001, Aharonian 2000). The relative contributions of 
various classes of AGN to the neutrino, gamma-ray and cosmic ray backgrounds are 
discussed by Mannheim et al. (2001). 

Interestingly, Tinyakov & Tkachev (2001b, 2003) have claimed a correlation between 
the arrival directions of cosmic rays, having energies above 240 EeV from the Yakutsk 
array and above 480 EeV from the AGASA array, with the the directions of the 22 most 
powerful BL Lac objects with redshifts z>Q.l. In this analysis, the energy cuts were 
those for which there was an indication of small-angle clustering from their previous 
analysis (Tinyakov & Tkachev 2001a). As pointed out by Evans, Ferrer & Sarkar (2003) 
the redshift cut implies that cosmic rays from these sources would be strongly affected 
by the GZK cut-off. Furthermore, cosmic rays propagating to Earth from large redshifts 
through the extragalactic magnetic field and CMBR may have difficulty in reaching us 
in a Hubble time as well as being severely affected by photoproduction interactions. 
Evans, Ferrer & Sarkar (2003) fail to find a statistically significant correlation between 
the cosmic ray arrival directions and the directions of BL Lac objects, and claim the 
correlation found by Tinyakov & Tkachev (2001b, 2003) to be spurious, and due to the 
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cuts imposed on the data. We would add that there is no reason to hmit searches to 
BL Lac objects as these objects are thought to be Fanaroff Riley Class I radio galaxies 
having jets closely ahgned to our hne of sight - unless the particles being searched for 
are neutral particles produced and relativistically beamed in the jets, there is no reason 
to favour BL Lac objects over the more numerous Fanaroff Riley Class 1 radio galaxies. 
The redshift cut (2;>0.1) was used to ensure the BL Lac objects were powerful, but 
more powerful Fanaroff Riley Class I radio galaxies can appear less luminous than less 
powerful BL Lac objects, due to the absence of strong relativistic beaming of photons 
toward the observer. Excluding nearer BL Lac objects and Fanaroff Riley Class I radio 
galaxies seems to us to be illogical. Indeed, the sub-parsec scale jets of the nearby 
Fanaroff-Riley Class I radio galaxy M87, near the centre of the Virgo Cluster, may 
be a mis-aligned blazar of the BL Lac type, and may accelerate protons to ultra-high 
energies which could propagate to Earth if the magnetic field topology between our 
Galaxy and the Virgo Cluster were favourable (Protheroe, Donea and Reimer 2003). 

Quasars are also an interesting possibility. Farrar and Biermann (1998) found a 
correlation between the arrival directions of the five highest energy cosmic rays hav- 
ing well measured arrival directions and radio-loud flat spectrum radio quasars with 
redshifts ranging from 0.3 to 2.2. The probability of obtaining the observed correla- 
tion by chance is estimated to be 0.5%. Although the statistical significance is not 
overwhelming, and indeed other researchers find no statistically significant evidence 
for such a correlation (Hoffman 1999, Sigl et al. 2001), if further evidence is provided 
for this correlation the consequences would be far reaching. The distances to these 
AGN are far in excess of the energy-loss distance for pion photoproduction by protons. 
Furthermore, given the existence of intergalactic magnetic fields, any charged particle 
would be significantly defiected and there should be no arrival direction correlation 
with objects at such distances. Hence, the particles responsible would need to be be 
stable, neutral, and have a very low cross section for interaction with radiation. Of 
currently known particles, only neutrinos fit this description, however supersymmetric 
particles are another possibility, as suggested by Farrar and Biermann (1998), although 
this possibility now appears to be ruled out (Gorbunov, Raffelt & Semikoz 2001). 

7.1 Z-Bursts 

The difficulty of having high energy neutrinos producing the highest energy cosmic 
rays directly is circumvented if the neutrinos interact well before reaching Earth and 
produce a particle or particles which will produce a normal looking air shower. As 
suggested by Weiler (1999), this may occur due to interactions with the 1.9 K cosmic 
background neutrinos (see also Gelmini and Kusenko 2000). The clustering of rehc 
neutrinos in hot dark matter galactic halos would give an even denser nearby target 
for ultra high energy neutrinos as suggested by Fargion et al. (1999). The cross section 
is much larger for resonant production which would occur for a UHE neutrino with 
energy ~ 111^/2111^ ~ 4 x 10^^/ (m^/eV) eV. From the recent SNO results( Ahmad et al., 
2002), the Super-Kamiokande atmospheric neutrino results (Toshito et al., 2001) and 
the tritium P decay results (Bonn et al., 2001), Ahmad et al. (2002) concluded that the 
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sum of mass eigenvalues of oscillating neutrinos was in the range 0.05-8.4 eV. Hence, 
the required UHE neutrino energy for resonant Z° production is very interestingly at 
or above the energies of UHE CR near the GZK cut-off. In this "Z-Burst" scenario, 
the Z'^ produced would have a comparable energy to the UHE neutrino, and decay 
into leptons and hadrons (including nucleons) which could be detected as UHE CR. In 
principle, there is also the possibility of the determination of absolute neutrino masses 
from Z-bursts if the galactic cosmic ray spectrum from normal acceleration were known 
(Pas and Weiler 2001, Fodor et al. 2002). 

The main problem with the Z-Burst scenarios is that except for the case of unreal- 
istic source models (production of UHE neutrinos with very few other UHE particles) 
or rather extreme over-densities (by >10^) of relic neutrinos, the cascade gamma-ray 
flux would exceed the GeV gamma-ray intensity observed by EGRET (e.g. Kalashev 
et al. 2002). Also, extremely high fluxes of UHE neutrinos are required to explain the 
observed UHE CR spectrum. They are well in excess of those expected from AGN, 
and the possibility of explanation in terms of X particle decay exclusively into neu- 
trinos seems unsatisfactory (Berezinsky et al. 2002). Nevertheless, these fluxes are in 
principle detectable with existing neutrino detectors, and if they exist should certainly 
be detected with future large area cosmic ray detectors, such as the Pierre Auger Ob- 
servatory, which are also sensitive to UHE neutrinos. If detected, the observed UHE 
neutrino flux, together with the observed UHE CR flux and anisotropy would place 
important constraints on the mass of the possible neutrino species forming the dark 
matter galactic halo, and its radial distribution (e.g. Singh and Ma 2003) as well as 
the source of the UHE neutrinos. Their detection would also require a re-evaluation of 
the our understanding of electromagnetic and hadronic cascading in the CMBR and 
other radiation flelds. 

8 Topological Defects 

Topological defects (TD) such as monopoles, cosmic strings, monopoles connected by 
strings, etc., may be produced at the post-inflation stage of the early Universe. In 
the process of their evolution the constituent superheavy flelds (particles) may be 
emitted through cusps of superconducting strings, during annihilation of monopole- 
antimonopolc pairs, etc. These particles, collectively called X-particles, can be su- 
perheavy Higgs particles, gauge bosons and massive supersymmetric (SUSY) parti- 
cles. These are generally very short-hved, and their decay followed by a hadroniza- 
tion cascade could produce an observable signal. Signals of TD origin would be af- 
fected by interactions/cascading during propagation over cosmological distances to 
Earth. Protheroe and Johnson (1996) pointed out the importance of including pair- 
synchrotron cascades in UHE CR propagation and, following their approach, Protheroe 
and Stanev (1996) showed that the 7-ray flux for many TD models of UHE CR exceeded 
that observed at 100 MeV energies for B > 10'^ G. 

There could be also superheavy quasi-stable particles with lifetimes larger (or much 
larger) than the age of the Universe. These particles could be produced by many 
mechanisms during the post-inflation epoch, and survive until the present epoch. One 



32 



interesting process is the "gravitational production of super-heavy particles" , in which 
no interaction of X-particles is required. Also, string theories predict the existence of 
other super-heavy particles ("cryptons") which are metastable and could in principle 
form part of the cold dark matter (CDM) (see, e.g., Kolb 1998, Ellis 2000). As with 
any other kind of CDM, super- heavy quasi-stable X-particles would cluster in galactic 
halos. The same clustering would also occur for some TD, such as monopolonium, 
monopole-antimonopole pairs connected by a string, and vortons. Cosmic ray signals 
from all these objects would reach us relatively attenuated. Perhaps the most promising 
WIMP CDM candidate is the lightest SUSY particle (LSP) with mass only 20-1000 
GeV, and so would not produce UHE CR. 



8.1 Fragmentation functions 



TD such as cosmic strings, necklaces, etc., are extragalactic, and could produce an ex- 
tragalactic signal through the decay of short-lived X-particles. TD which accumulate in 
galaxy halos (monopolonia, monopole-antimonopole-pairs and vortons) could produce 
a galactic signal through the annihilation/emission and decay of short lived X-particles 
which would in turn decay promptly into Standard Model (SM) states. Super-heavy 
quasi-stable particles (r ^ to) would decay similarly, and also be clustered as CDM in 
galactic halos. 

The X particle decay products annihilate and could give rise to a jet of hadrons. 



e.g. 
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Energy spectra of the emerging particles, the "fragmentation functions", were first 
calculated by Hill (1983). Each jet has energy mxc^/2, and so one defines a dimen- 
sionless energy for the cascade particles, x = 2E /mxc^- The fragmentation function 
for "species a" is then defined as dNa/dx. A very fiat spectrum of particles results, and 
extends up to ~ mxc^/'2. In the case of decay of CDM in galactic halos, the resulting 
UHE CR spectrum is proportional the fragmentation function for nucleons. 

Some recent calculations of the fragmentation functions used the Modified Leading 
Logarithm Approximation (MLLA) which is valid only for a; <^ 1, and in more recent 
QCD calculations PYTHIA/JETSET (Singh and Ma 2003) or HERWIG Monte Carlo 
event generators were used. The fragmentation functions due to Hill (1983), and those 
of Berezinsky et al. (1997) based on the MLLA are compared in Fig. 14. Initially, 
the inclusion of the production of SUSY particles was done by putting 40% of the 
cascade energy above threshold for production of SUSY particles into LSP (Berezin- 
sky and Kachelricss 1998), thereby steepening the fragmentation functions for normal 
particles at high energy. Birkel and Sarkar (1998) showed that even without inclusion 
of SUSY production there is a significant dependence on nix, such that for high mx 
the fragmentation functions are steeper, as a direct consequence of the well-known 
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Figure 14: Fragmentation functions for hadronization of baryons normalized at x = 0.1 
to the recent QCD calculation of Sarkar and Toldra (2002) for nix = lO^^GeV (solid 
curve): Hill(1983) (dot-dashed curve); Berezinsky et al. (1997) MLLA approxima- 
tion (short dashed curve); Monte Carlo results of Birkel and Sarkar (1998) for mx ~ 
lO^^GcV (solid histogram); Berezinsky and Kachclriess (2000) for mx = lO^^GeV 
SUSY-QCD (long dashed curve); Sarkar and Toldra (2002) for mx = lO^^GeV SUSY- 
QCD (dotted curve); Rubin (2000) mx = IO^^q^v QCD (dot-dot-dot-dashed curve). 

Feynman scaling violation in QCD. Fragmentation functions calculated by Birkel and 
Sarkar (1998) using the HERWIG event generator have been added to Fig. 14, but this 

event generator is now known to overestimate production of nucleons by a factor ~2- 
3, (Sarkar 2000, Rubin 2000). Recent calculations by Rubin (2000), Sarkar and Toldra 
(2002) and Berezinsky and Kachelriess (2000) (added to Fig. 14) have used improved 
treatments of SUSY particle production, and result in only ~ 5-12% of the cascade 
energy going into LSP. Very recently, Barbot & Drees (2003) have produced a complete 
set of fragmentation functions for any SUSY particle of the minimal supersymmetric 
extension of the standard model into protons, photons, electrons, neutrinos and the 
LSR 

8.2 Viability of Dark Matter Origin of UHECR 

Predictions for some TD and massive relic particles are show in Fig. 15. If CDM consists 
of particles associated with TD distributed uniformly throughout the Universe, then 
UHE CR are subject to the GZK cut-off. In this case 7-ray signals result from a pair- 
synchrotron cascade in background radiation and extragalactic magnetic fields. The 
magnetic fields used in some cascade calculations may have been unrealistically low. 
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Figure 15: Predicted TD model spectra of protons (thick curves) and cascade photons 
(thin curves) compared with cosmic ray data from Gaisser and Stanev (1998). (a) 
mx = 10^^ GeV: Protheroe & Stanev (1996) (sohd curves); Sigl et ah (1999) X ^ uu 
(short dashed curves); Blasi (1999) super-heavy rehc halo population, SUSY-QCD and 
TT/xe synchrotron (gamma rays only, long dashed curve), QCD and vr/xe synchrotron 
(gamma rays only, dot-dot-dot-dashed curve); Berezinsky et al. (1998) super-heavy 
relic halo population, SUSY-QCD (dotted curves), necklaces SUSY-QCD (dot-dashed 
curves), (b) mx = lO^^-W^ GeV: Sigl et al. (1999) mx = lO^^ Qgy, x ^ q + q 
QCD (solid curves), mx = 10^^ Qgy, X ^ q + l QCD (dotted curves); Blasi (1999) 
mx = 10^^ GeV super-heavy relic halo population, SUSY-QCD and iifj,e synchrotron 
(gamma rays only, long dashed curve), QCD and vr/ie synchrotron (gamma rays only, 
dot-dot-dot-dashed curve); Sarkar and Toldra (2002) super- heavy rehc halo population 
QCD best fit mx = 10^^ GeV (protons only, leftmost dot-dashed curve) and SUSY- 
QCD best fit mx = 5 x 10^^ GeV (rightmost dot-dashed curve). 

and it does appear difficult to explain the super-GZK events with such TD models 
without the flux of cascade gamma-rays exceeding the observed 100 MeV gamma-ray 
background. 

Most of the matter in the Universe is CDM, and if it consists of massive relic 
particles they should cluster in galaxy halos. In this case, decay of massive relic particles 
would produce UHECR signals weakly anisotropic toward the GC, and the UHE CR 
spectrum would not have a GZK cut-off. Of the "top down" scenarios, these models 
currently seem to us to be the most viable. See Sarkar and Toldra (2002) for recent 
work on decay of superheavy dark matter particles. 
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9 Propagation through Magnetic Fields 



Cosmic rays reach us after travelling through the magnetic fields which pervade space. 
Details of the strength and structure of such fields are unknown but broad generaliza- 
tions are possible within certain volumes of the Universe. 

Our galaxy is of spiral structure and the galactic magnetic field has a regular com- 
ponent with a characteristic strength of the order of microgauss and which seems to be 
associated with the spiral arms. Additionally, there is a turbulent, random, component 
which is at least as strong as the regular component. This component is even less well 
known since measurements of Faraday rotation or other techniques tend to average 
out when the line of sight transits a number of turbulence cells. The spectrum of tur- 
bulence scales is often assumed to be of a Kolmogorov kind which has the important 
property of being dominated by the largest scale sizes. Within our Galaxy, the largest 
internal structures tend to be of 100 pc scales (e.g. supernova remnants). With a field 
strength of a few microgauss, this means that significant scattering of cosmic rays will 
occur at least to a few times 10^'' eV since a proton with this energy has a gyroradius 
of 100 pc in a 1 //G field. Hence, the largest scale lengths in the turbulence of the 
Galactic magnetic field tend to dominate the propagation of the (more energetic) UHE 
CR in the Galaxy. 

The particles follow paths rather like random walks up to the scale size of the 
turbulence. A first approximation to galactic propagation is then diffusion. Honda 
(1987) has given an excellent discussion of extensions to make the simple picture more 
realistic. That work, and similar propagation modelling by Clay (2000), gives us some 
understanding of resulting measurable properties of the cosmic ray beam. Since the 
propagation is diffusive, the time for a particle to leave the galaxy is greater than 
the simple direct transit time. This containment time increase results in an increase 
in flux over that which would have been observed had there only been straight line 
propagation from galactic sources. Containment time calculations thus allow us to 
crudely determine a 'source spectrum'. In one recent analysis (Clay 2000), the source 
spectrum shows no knee and, possibly, no ankle. It may be that both those features 
(and certainly the knee) are consistent with purely propagation effects. The resulting 
source spectrum is a power law with an index of 2. There is a problem with such an 
explanation for the ankle since particles above that energy travel in rather straight 
lines and a 'Milky Way' perhaps ought to be visible in the anisotropy data. However, 
data at these energies arc somewhat sparse and, of course, the AGAS A/SUGAR source 
(Hayashida et al. 1999, Bellido et al. 2001) could be just such an effect. 

In considering the propagation from extragalactic sources, there are two environ- 
ments to consider. They are the intra-cluster magnetic flelds of both the source galaxy 
and our own galaxy, and the inter-cluster field. Clarke ct al. (2001) have shown that, 
remarkably, a characteristic intra-cluster magnetic field strength in a rich galactic clus- 
ter fills the cluster and has microgauss strengths - maybe 5/iG in the inner 500 kpc. If 
we make a first approximation to a diffusion coefficient to be a factor 77 > 1 times the 
minimum (Bohm) diffusion coefficient, i.e. r)rgc/3, we can assume diffusive propagation 
and derive an estimate of the time to reach a given root- mean-square displacement. If 
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Figure 16: Angular distribution of the arrival angle at Earth for mono-energetic injec- 
tion of protons of energy E ~ 10^^'^ eV, and for various source distances as indicated 
(from Stanev et al. 2000). 

we consider a lO^^cV particle in the 5/xG field, we find a required time of 10^/?] yr just 
to leave the source cluster through 500 kpc. The GZK effect is clearly relevant here. 
If the source is in the Virgo Cluster of galaxies, a cosmic ray must then travel through 
intercluster space and then reach us though our own cluster field. It may be that the 
intcrcluster field is also at significant levels. In this case, we are looking at tens of 
mcgaparsec from the nearest likely AGN source with a transit time, being dependent 
on the square of the distance, and becoming greater than the age of the Universe. This 
is clearly an issue which pushes us to a careful consideration of very local sources. 

This argument is rather crude. Sigl (2000) has modelled time delays for particles 
travelling 10 Mpc in a turbulent 0.3/iG field. Even for that modest field strength, 
transit times of 10^ yr apply at 10-100 EeV. Apart from any concern about particles 
reaching us within the age of the Universe, our comments on diffusion times emphasize 
that it may not make sense to correlate cosmic ray observations with sources beyond 
10 Mpc unless one can be sure that those sources have a lifetime for emission sub- 
stantially greater than lO^yr. Monte Carlo calculations have been made by Stanev et 
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al. (2000) for propagation through an irregular magnetic field having a Kolmogorov 
spectrum of turbulence with minimum wavenumber 1 Mpc~^ and energy density equal 
to that of a much smaller 1 nG intergalactic field. Propagation from sources at several 
distances take account of diffusion in this turbulent field as well as interactions with 
the CMBR using the SOPHIA event generator (Miicke et al. 2000) and rcdshifting, and 
give the distributions in energy, time delay relative to straight-line propagation, and 
angular distribution about source direction. Fig. 16 shows the angular distribution of 
protons with initial energies 300 EeV arriving from sources at distances 2 Mpc, 8 Mpc, 
. . . 512 Mpc away. 

The situation is quite different if the intergalactic magnetic field structure is based 
on the observation of microgauss fields in clusters of galaxies (Kronberg 1994), and of 
clusters occurring in networks of "walls" separated by "voids" . Using a wall/void model 
similar to that of Medina Tanco (1998), Protheroe et al. (2003) discuss propagation 
from M87 (located close to the centre of the Virgo Cluster) assuming that it and our 
galaxy are embedded in the same wall (thickness 2.5 Mpc) which has a regular magnetic 
field of 10~^ G in the plane of the wall and 10~^° G in the surrounding void, and an 
irregular component with 30% of the energy density of the regular component and 
having a Kolmogorov spectrum of turbulence. ModeUing M87 as a mis-aligned BL Lac 
object, they found the UHECR output from M87 to be at a level such that if UHECRs 
travelled in straight lines they would give an average intensity at Earth a factor ~20 
below that observed. However, observed fiuxes will be very different for propagation 
in a more realistic magnetic field structure such as that described above. Propagation 
results for the wall/void model are shown for two initial energies 10^^ eV and 10^° eV 
in Fig. 17. Note that Fig. 17 does not give the anisotropy that would be observed at 
the Earth due to cosmic rays from M87, but rather the enhancement factor (relative to 
straight line propagation) of the cosmic ray fiux from M87 at positions on a sphere of 
radius 16 Mpc centred on M87. Enhancement factors ~10^ exist if the Earth is within 
~1.5 Mpc of a field fine originating at M87. If this (rather special) condition were 
met M87 could easily explain the observed UHECR. Because of its very high black 
hole mass M87 was probably much more active at earlier times than at present (many 
objects exhibit a high state for ~5% of the time) which makes M87 a more attractive 
candidate source of the UHECR. 

10 Conclusions 

Cosmic ray astrophysics at the highest energies is entering a new era with major new 
facilities being proposed and developed. The tiny flux of particles in this range of 
energies challenges our understanding of the astrophysics which we apply to many 
other studies. We are led to examine processes in extreme regions such as the various 
components of AGN, to ask how particle acceleration may take place there, and to ask 
for detailed information on photon and magnetic fields throughout our local universe. 
This challenging work now truly complements studies throughout photon astrophysics. 
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Figure 17: Enhancement factor relative to rectilinear propagation as a function position 
(cos 9,(1)) a sphere of radius ~ 16 Mpc centred on the origin for protons of initial energy 
(a) Eq = 10^^ eV and (b) Eq = 10^° eV injected isotropically at the origin in the 
wall/void model of the IGMF discussed in the text. Positions (0,0°) and (0, 180°) are 
on the field line threading M87 (from Protheroe et al. 2003). 
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